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Abstract 

We explore the possibility of using the method of classical integral transforms to 
solve a class of g-difference-differential equations. The Laplace and the Mellin transform 
of (/-derivatives are derived. The results show that the Mellin transform of the q- 
derivative resembles most closely the corresponding expression in classical analysis, 
and it could therefore be useful in solving certain ^-difference equations. 



Revised version 



1 



1. The study of g- analysis is an old subject, which dates back to the end of the 19th 
century ([Q-Q]). It has found many applications in such areas as the theory of partitions, 
combinatorics, exactly solvable models in statistical mechanics, computer algebra, etc 
Recent developments in the theory of quantum group has boosted further interests in this 
old subject [|], 0. 

The subject of g-analysis concerns mainly the properties of the so-called g-special func- 
tions, which are the extensions of the classical special functions based on a parameter, or 
the base, q. The relations among these functions, and the difference equations satisfied by 
them are among the topics most studied so far. The g-difference equations involve a new 
kind of difference operator, the g-derivative, which can be viewed as a sort of deformation 
of the ordinary derivative. Solutions of the g-difference equations in one variable have been 
well studied in terms of the g-hypergeometric series (also called the basic hypergeometric 
series). Partial g-difference equations and g-difference-differential equations with more than 
one variables are generally studied by means of the method of seperation of variables, or 
by the techniques of Lie symmetry in the literature ([0,@|-[0I)- '^^^ method of integral 
transforms, which is another powerful technique of solving differential equations in classical 
analysis, has not been, in our view, explored in g-analysis. The reason is not hard to under- 
stand. The main virtue of the classical integral transforms, particularly the Fourier and the 
Laplace transform, is to transform a differential equation into an algebraic equation, which 
can be solved easily. That this is possible is due to the fact that these transforms change 
the derivatives of a function to something proportional to the transform of the original func- 
tion. As far as we know, integral tranforms or g-integral transforms which could transform 
g-difference equations into algebraic equations have not been found. It should be mentioned 
that in fact g-analogues of Fourier transform, based on the Jackson g-integral, have been 
proposed recently [l^, |15||. However, in order for the g- Fourier transform of the g-derivative 
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of a function f{x) to be proportional to the g- Fourier transform of f{x), the function f{x) 
must satisfy very special conditions, such as f{q^^) = = f{—q^^) Ijl^- Hence, while 



these g-Fourier transforms may be useful in proving certain identities among the g-special 
functions, their use in solving g-difference equations seems limited. 

In this paper we shall explore the possibility of using the method of classical integral 
transform to solve a class of g-difference-differential equations. We derive the Laplace and 
the Mellin transform of g-derivative, and argue that the Mellin transform, which is not 
generally employed in solving differential equations in classical analysis, may still be useful 
in solving certain g-difference equations. 

2. Suppose we want to solve the following g-diffusion equation 

D^y{x,t) = —y{x,t) (-00 < x < 00, t > 0) (1) 
subject to the initial condition 

y{x,0) = f{x). (2) 



Here Df is the "forward" temporal g-derivative defined by ^ 



ntm '^^j^ ■ (3) 

for any function h{x). We assume < g < 1 in this paper. The function f{x) is assumed to 
vanish as x — ±00. One may as well use the more common definition of g-derivative 0] 

VMt) := . (4) 

(1 -g)t 

We shall not employ this definition of the g-derivative here for reason to be explained later. 
We note here that g-difference and g-difference-differential equations of the diffusion type 
such as eq.(|l|) have been considered before but mostly from the point of view of 

Lie symmetry, or by seperation of variables. 



We can remove the partial differential operator in x in (0) by a Fourier transform. The 
question now is to choose an appropriate integral transform to remove the g-derivative. In 
view of the positivity of the time variable, the two most natural choices are the Laplace and 
the Mellin transform. 

Let us first derive the expression of the Laplace transform of the g-derivative. The Laplace 
transform of a function h(t) is defined as h{s) := C{h{x),s} = h{t) exp{—st)dt. For the 
g-derivative of h{x), the Laplace transform is 

1 



C{D!hit),s} 



t 



t 



(5) 



To proceed we have to use the following relation of the Laplace transform [p!S 

h{s')ds' 



hit) 



e-''dt , 



(6) 
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provided the integral on the r.h.s. of (||) is well-defined. We may apply ® to (|]) directly if 
/i(0) = 0. However, if h{{]) ^ 0, the r.h.s. of (P) is not well-defined, and direct application 
of (I) to @ leads to incorrect result which does not reduce to the usual expression of the 
Laplace transform of derivative in the classical limit g — 1~. In order to recover the classical 
limit correctly, we find it necessary to regularise (|]) in the form 
1 



1-g 



%-H)-M0)^-.,,_ 



(7) 



We may now apply (D to (0). Making use of 



C{h{t) - /i(0), s} = h{s) - s-^h{0) 



we finally obtained 



C{Dlh{t),s} = ^ r Hs) ds' - ^/i(O) . 

i — Q J sq i — g 



(9) 



Eq.(P) reduces to the expression sh{s) — h{0) for the Laplace transform of ordinary derivative 



If one uses instead the definition @ for the g-derivative, the Laplace transform would be 

CiVlKt), s} = -^ h{s') ds' - ^-^h{0) . (10) 
I — q Js 1 — ? 

It is now obvious that the Laplace transform is not useful in solving equations involving 

g-derivatives: it transforms such equations into integral equations! 

3. We now consider the Mellin transform of a g-derivative. The Mellin transform is 
seldom being used in solving differential equations, because it generally transforms differ- 
ential equations into difference equations instead of the much simpler algebraic equations. 
Now that the Fourier and the Laplace transform lose their virtues whenever g-derivatives 
are present, the Mellin transform is naturally the next one to be looked at. As we shall 
see below, the Mellin transform still transforms an equation containing g-derivatives into a 
difference equation of the transformed function, which is the best thing next to an algebraic 
equation one could get. Previously, the use of the Mellin transform in g-analysis is limited 
to proving various identities among the g-special functions [|, 0. 

The Mellin transform of a function h(t) is defined as h*{s) := M{h(t), s} = /g°° h{t)t''^^dt. 
For g-derivative defined in (^, we have 

M{D!h{t), s} = -[s- l]g h*{s - 1) . (11) 

Here [x]g is the g-number defined by 

:= . (12) 

Note that [ 1 . Hence ([Tl|) reduces to the expression — (s — l)h*{s — 1) for 

the Mellin transform of the ordinary derivative as g 1^. Repeated use of ([Tl|) leads to 

MiiDirhit), s} = i-iris - lUs ~2]g---[s-n],h*is-n) , n>l. (13) 



This is the g-analogue of the corresponding formula in the classical case 



For the definition one has 



M{Djh{t),s} = [l-.s],h*{s~l) (14) 
= -q^-^[s-l]^h*{s-l) . (15) 



Here an extra factor of q appears compared with (11). In order to simphfy our presentation, 
we therefore adopt the definition (§) in this paper. We must, however, mention that all 
the arguments given below apply equally well to the corresponding cases with g-derivatives 
replaced by the definition (|^). 

4. Let Y*{^,s) be the transformed function of y{x,t) obtained by taking the Mellin 
transform in t and a Fourier transform G{^) := f!^^ g{x) exp{iC,x)dx in x. Making these 
transforms to ([l|), one obtains 

[s-l],Y*i^,s-l) = eY*i^,s) . (16) 

Fortunately solution to this equation can be readily found to be 

Y*{^,s) = A{Or''T,{s) , (17) 

where A{^) is some function of C, only, and Tg{s) is the q-gamma function defined by 



(a;9)oo:= 11(1 -<"'')■ (19) 

k=0 

Tq{s) satisfies 

lim T,{s) = T{s) , (20) 
T,is + 1) = [s],Tis), r,(l) = l. (21) 
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Inverse- Mellin transform of ^ ^^rg(s) in (0) is 

TT-. / i-''T,{s)t-^ds . (22) 

The poles of Tq{s) are s = 0, —1, —2, . . .. The residual of Tq{s) at pole s = —n (n > 0) is [Q: 

(1 - 



(g-";g)„lng-i 
The symbol (a; g)„ is the q-shifted factorial: 



Hence (^) becomes 



In view of the identity [Q 



(23) 



(a;g)o:=l, n = 0, (24) 

{a-q)n := (1 - a) (1 - ag) ■ ■ ■ (1 - ag""^) , n=l,2... (25) 



(r^g)n= f--| g~'^("-')/'(g;g)n, (27) 



can be expressed as 

1-q ^ _ g)^2^]» 

= ^,Eq{-qil-q)et) . (28) 
The function Eq{z) (for complex 2;) is the q-exponential function defined by 

00 n{n—l)/2^n 

E,{z) := Y: \ . , = (-^; <l)oo ■ (29) 

In the limit g 1^, eq.(^) tends to the usual exponential function exp(— ,^^t). Finally, 
performing an inverse Fourier transform we obtain the solution of the g-diffusion equation 

yix, = 4= ^(0 I ^ Eq (-g(l - q)ex) \ e^^^ . (30) 



Setting t = in (0) shows that 



1-q 



In q~ 



km 



1 /"oo 



'271 J- 



y{x, 0)e*^^ dx 
f{x)e'^'^ dx 



(31) 



is the Fourier transform of y{x, 0) = f{x). So the final solution of the initial problem is 

y{x, t) = ^ FiOE, (-g(l - q)et) e'^^^ d^ . (32) 

This is the g-analogue of the solution given in for the corresponding classical case. One 
can easily check that (|3^) indeed satisfies (|l|) by using the following identity 



DlE^{\t) 



A 



-E,{\t) . 



(33) 



q{l - q) 

Let us consider an example. Suppose the initial profile is f{x) = exp(— a;^/46)/-\/26 
(6 > 0). Its Fourier transform is F{^) = exp(— 6^^). Then from (P2|) and (P^), we get 



y{x,t)=EJq{l-q)t-^]f{x 



(34) 



In the limit q ^ 1 , eq. (|3^) gives the classical solution 



y{x,t) 



e db 



1 _^ 

e 46 



. 26 

1 ^2 

g 4(t+6) 



V2(t + &) 

5. As another example, let us consider the following wave equation 



(35) 



{Dlfy{x,t) 



dx^ 



y{x,t) (— oo < X < oo , t > 0) 



(36) 



with inital conditions 



yix, 0) = /(x), A^?/(x, 0) = g{x) 



(37) 



We assume that both f{x) and g{x) vanish as x — > ±00. In this case the Fourier- Melhn 
transformed function Y*{^,s) obeys 

[s - lUs - 2],r*(e, s-2) = -eY*i^, s) . (38) 

The general solution is 

s) = [a{0 HO'' + BiO m^'] r,(5) , (39) 

where A{C,) and B{C,) are some functions of ^. Performing the inverse- Mellin transform, we 
get 

Y{^, t) = l^,{MOEg (^g(l - qm + B{i)E, {-iq{l - q)it)} . (40) 

Here Y{C,,t) is the Fourier transform of y{x,t) with respect to x. Now we rewrite (|40| ) in 
terms of the q-Sine and the q-Cosine function which are defined by @] 

Sin,(x) = ^^^^"^ -f ^^-^"^ , (41) 

Cos,(x) = ^^^^"^ . (42) 

The result is 

y{^, t) = ]^{C(0 Cos, (g(l - g)et) + D{0 Sin, (g(l - q)^t)] , (43) 

where the functions C(^) and -D(0 c'-^'g linear combinations of A{0 and -8(0- The inverse- 
Fourier transform of (|43|) is 

|/(x, t) = ^ip— {C(0 Cos, (g(l - qm + D{0 Sin, (g(l - q)^t)} e"^?^' . (44) 
V^TT In q~^ J -00 ^ 

Letting t = in (0), one can check that the function C(^) is related to the Fourier transform 
of f\x) by 

m = m . (45) 
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Making use of the following relations, which can be obtained by means of (|33| ) 



Dt Sin,(At) 



A 



COSg(At) 



g(l - q) 

Dl Cos,(At) = -^^X^ ^iri^iXt) , 
we can relate D{^) to the Fourier transform G{^) of g{x) as follows: 



(46) 
(47) 



G{i) = ^ Dm ■ 

With these results, we finally obtain the solution to the initial problem of eq.( 



(4J 



y{x,t) 



^27r J- 



F{0 Cos, (g(l - qm + ^ Sin, (g(l - q^] e''^^ 



(49) 



This solution is the g-analogue of the solution to the corresponding classical case given in 



19 . 



6. We now see how the above steps are generalised to the equation: 

{Dlf y{x,t) = — y{x,t) (-cx)< a; < oo , t > 0, n>2) 



(50) 



with inital conditions 



l/(a;,0) = /(x), {Dif y{x,Q) = gk{x) , = 1, . . . , n - 1 , 



(51) 



where the functions f{x) and gk{x) are assumed to vanish as x ^ ±oo. The Fourier-Mellin 
transformed function 1^*(^, s) obeys 



'-ins - lUs - 2], ■ ■ ■ [s - n],y*(e, s-n) = -^Y^i, s 



(52) 



The general solution is 



n-l 



F*(e,a) = r,(.)r^ ^ A^{i) 



m=0 



— e 



_(2m + l). 



(53) 
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where >lm(0 ^-^e some functions of ^. We can now perform the inverse MeUin and Fourier 
transforms to get the final solution, which is given formally as 

V27r In g-i J-oo V / 

The functions ^^(0 can then be related to the Fourier transforms of the functions f{x) and 
gk{x) from the initial conditions. 

7. To summarise, we show that the Mellin transform of the g-derivative resembles most 
closely the corresponding expression in classical analysis, whereas transforms such as the 
Fourier and the Laplace transform fail in this respect. As such the Mellin transform can be 
useful in solving certain g-difference equations. We illustrated this fact with a few examples. 
However, for the MeUin transform to be really useful, a more complete knowledge of the 
properties of the g-special functions under various integral transforms (Fourier, Laplace, 
MeUin, etc) and their inverses has yet to be attained. What is more desirable is to invent 
integral transforms or g-integral transforms that possess the virtue of the Fourier and the 
Laplace transform in the classical analysis mentioned in the introduction. 
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